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Abstract. We explore the connection between the transfer matrix formalism and discrete complex 
analysis approach to the two dimensional Ising model. 

We construct a discrete analytic continuation matrix, analyze its spectrum and establish a direct 
connection with the critical Ising transfer matrix. We show that the lattice fermion operators of 
the transfer matrix formalism satisfy, as operators, discrete holomorphicity, and we show that their 
correlation functions are Ising parafermionic observables. We extend these correspondences also to 
outside the critical point. 

We show that critical Ising correlations can be computed with operators on discrete Cauchy data 
spaces, which encode the geometry and operator insertions in a manner analogous to the quantum 
states in the transfer matrix formalism. 



(N 



> 
o 

o 

CO 
Oh 

i Contents 

1. Introduction 2 

1.1. Ising model 2 
i_i 1.2. Transfer matrix and discrete holomorphicity 3 

1.3. Main results 4 

2. S-Holomorphicity and Riemann Boundary values 7 
Q\ 2.1. S-holomorphicity equations 7 

2.2. Massive s-holomorphicity 8 

2.3. Riemann boundary values 9 

2.4. S-holomorphic continuation operator 9 

2.5. Spectral splitting of the propagator 11 
fN| 3. Transfer Matrix, Clifford Algebra and Induced Rotation 12 

3.1. Transfer matrix and Clifford algebra 12 

>• 3.2. Induced rotation 13 

3.3. Fock representations 15 

3.4. A simple polarization for low temperature expansions 16 

3.5. The physical polarization 17 

4. Operator Correlations and Observables 18 

4.1. Operator insertions in the Ising model transfer matrix formalism 18 

4.2. S-holomorphicity of fermion operator 20 

4.3. Ising parafermionic observables and low temperature expansions 21 

4.4. Pfafnan formulas for multi-point fermion correlation functions 24 

4.5. Multipoint Ising parafermionic observables 25 

4.6. Correlation functions of the fermion and spin operators 26 

5. Operators on Cauchy data spaces 27 

5.1. Discrete RPS operators 27 

5.2. RPS pairings 28 

5.3. Fermion correlation and fixed point problems Cauchy data 29 
References 30 



DISCRETE HOLOMORPHICITY AND ISING MODEL OPERATOR FORMALISM 



2 



1. Introduction 

The transfer matrix approach to the planar Ising model is both classical and remarkably powerful 
[KrWa41, Bax82, McWu73]. The free energy, critical exponents and a number of correlation func- 
tions of the model were calculated using the transfer matrix, and much of the algebraic structure 
underlying the model is easiest understood by means of the transfer matrix and the related operator 
formalism. The formalism is also manifestly suggestive of the quantum field theories believed to 
describe the scaling limit of the Ising model. 

Recently, methods of discrete complex analysis have lead to significant progress in the under- 
standing of the Ising model, especially in its critical phase [Smi06, SmilOb]. The discrete complex 
analysis techniques apply to the model on planar domains of arbitrary shapes, and allow to prove 
conformal invariance results. 

In this paper, we investigate connections between the transfer matrix approach and discrete 
complex analysis techniques. We study discrete level relations between the two approaches, using 
concepts of quantum field theory and analytic tools that are well behaved in the scaling limit. 

1.1. Ising model. The Ising model describes up/down spins interacting on a lattice. It is a simple 
model originally introduced to describe ferromagnetism, but subsequently it has become a standard 
in the study of order-disorder phase transition. 

The Ising model is a random assignment of ±1 spins to the vertices a graph, that interact via the 
edges of the graph. We will consider the Ising model on subgraphs Q = (V, £ ) of the square lattice 
I?. 

The probability of a spin configuration (s x ) xgV is proportional to the Boltzmann weight e~@ H ( s \ 
where /3 > is the inverse temperature and H is the energy given by H(s) = — s i s j with sum 
over pairs of adjacent vertices < i,j >G £. Hence, the model favors local alignment of spins by 
assigning them a lower energy, and the strength of this effect is controlled by f3. 

For the Ising model in dimensions at least two, a phase transition in the large scale behavior 
occurs at a critical inverse temperature, for the square lattice Ising model at (5 C = |ln(\/2 + l). 
For (3 < f3 c , the system is disordered: spins at large distances decorrelate, i.e. there is no alignment. 
For (3 > (3 C the system has long-range order: spins are uniformly positively correlated, i.e. global 
alignment takes place. To properly make sense of the large scale behavior, one considers either the 
thermodynamic limit in which the graph tends to the infinite square lattice Q / Z 2 , or the scaling 
limit in which a given planar domain Q, is approximated by subgraphs £1$ °f SZ, 2 , the square lattice 
with fine lattice mesh 5 \ 0. 

The square-lattice Ising model is exactly solvable: in particular, the free energy and thermody- 
namical properties of the model are well understood. However, the fine nature of the critical phase 
and its precise connection with quantum field theory have for long remained mysterious from a 
mathematical perspective. Renormalization group and quantum field theory methods have provided 
a non-rigorous insight into the nature of the phase transition, Conformal Field Theory in particular 
giving numerous exact predictions. At the critical point, (3 = f3 c , the model (like many critical two 
dimensional lattice models) should have a universal, conformally invariant scaling limit. Recently 
some of this insight has become tractable mathematically by the development of discrete complex 
analysis techniques: one can make sense of the scaling limits of the fields and the curves of the model 
at the critical temperature. 

• The scaling limits of the random fields of the model are described by a Conformal Field 
Theory. The CFTs are quantum field theories with infinite-dimensional symmetries, which 
allow one to compute the critical exponents and the correlation functions via representation 
theoretic methods [BPZ84a, BPZ84b]. 

• The scaling limits of the random curves of the model are described by a Schramm-Loewner 
Evolution. The SLEs are random processes characterized by their conformal invariance and 
a Markovian property with respect to the domain [SchOO]. 
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A natural framework to investigate full conformal invariance of the Ising model (and other mod- 
els) is to study the model on arbitrary planar domains, with boundary conditions. A number 
of results in this framework has been obtained in recent years: the convergence in the scaling 
limit has been shown for parafermionic observables [Smi06, SmilOa, ChSm09], for the energy cor- 
relations [HoSmlOb, HonlOa] and for the spin correlations [Chlzll, CHI12]. These scaling limit 
results for correlations rely, for a large part, on discrete complex analysis. They have in turn pro- 
vided the key tools to identify and control convergence of the random curves in the scaling limit 
[Smi06, CDHKS12, HoKyll]. 

In the special case of the full plane, the progress in the study of scaling limits of Ising model at 
and near criticality has been steady over a longer time. In notable results, massive correlations in 
the full plane have been computed [WMTB76] and formulated in terms of holonomic field theory 
[SMJ77, SMJ79a, SMJ79b, SMJ80, PaTr83, Pal07]. Critical correlations in the full plane have been 
computed using dimer techniques and discrete analysis [BoDT09, B0DTO8, Dublla, Dubllb]. 

1.2. Transfer matrix and discrete holomorphicity. In this subsection, we briefly introduce 
the two approaches to the Ising model studied in this paper: the transfer matrix and the discrete 
complex analysis formalisms. 

1.2.1. Transfer matrix approach. Let I be an interval of Z, with boundary dl C I consisting of the 
two endpoints of the interval, and consider the rectangular box I x {0, . . . , N} with rows Io, . . . , I at 
(set l y := I X {y}). Using the transfer matrix, we can represent the Ising model on I x {0, ... , N} 
as a quantum evolution of spins living on I from time to time N. 

The Ising model transfer matrix V acts on a state space S, which has basis (e CT ) indexed by spin 

configurations in a row, a G {il} 1 - We set V := (V h )^ V v (U h ) 5 , where the factors separately 
account for Ising interactions along horizontal and vertical edges. The matrix element of V v at 
a, p G {il} 1 is defined (with fixed boundary conditions) as 



Viewing the y-axis as time, the transfer matrix can be thought of as an exponentiated quantum 
Hamiltonian in 1 + 1 dimensional space-time: at the row 1^, we have a quantum state v y G 5 which 
we propagate to the next row 1 y +i by v^+i = V\ y . 

In the path integral picture, the evolution vo, • • • , v^v becomes a sum over trajectories weighted 
by their amplitudes: the trajectories are spin configurations s : I x {0, ...,N} — > {±1} and the 
amplitudes are the Boltzmann weights e~@ H ( s \ As a result, the partition function Z = ^2 s e~^ H ^ 
equals (fjV^li) := f V^i, where i, f G S encode the boundary conditions on Io and Ijy. 

Ising fields (such as the spin, energy, disorder, fermions) are represented by the insertion of 
corresponding operators. The position of the fields appear in two ways: the operator is applied to 
the state on the row y on which the field lives, and the applied operator O x : S — > S depends on 
the position x of the field in that row. We combine the dependence on the horizontal coordinate x 
and the vertical time coordinate y by using the operator 0{z) = V~ y O x V y for the field located at 
z = x + \y. Then the correlation function of fields . . . , located at z±, . . . , z n is 




and the matrix V h 



is diagonal with elements 




(f\V N QW (zi) • • • OW ( Zn ) |i) 
(f\V N \i) 
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For probabilistic fields such as the spin s z , represented by a (z) : S — > S, the correlation functions 
are the expected values of products, e.g. (a (z\) ■ ■ ■ a (z n )) = E [s Zl • • • s z J. Non-probabilistic fields 
(such as fermion and disorder) can also be represented within the transfer matrix formalism. 

Also in Conformal Field Theory, the field-to-operator correspondence is fundamental. However, 
naively connecting the algebraic structure of Ising model and the one of CFT is problematic: the 
transfer matrix does not have a nice scaling limit and it is best suited to very specific geometries 
(rectangle, cylinder, torus, plane). 

Contrary to the transfer matrix formalism, discrete complex analysis is well suited to handle 
scaling limits on domains of arbitrary geometry, and hence to discuss conformal invariance. For this 
reason, relating the transfer matrix to discrete complex analysis seems a promising way to provide a 
manageable scaling limit for the quantum field theoretic concepts of the transfer matrix formalism. 

1.2.2. Discrete complex analysis approach. The idea of discrete complex analysis is to identify fields 
on lattice level, whose correlations satisfy difference equations — lattice analogues of equations of 
motion. A particularly useful type of such equations are strong lattice Cauchy-Riemann equations 
(massless at j3 c , massive at (3 / /3 C ), which we will refer to as ' s-holomorphicity 1 . 

For the critical Ising model, certain s-holomorphic fields can be completely characterized in terms 
of discrete complex analysis: their correlation functions (called 'observables') can be formulated as 
the unique solutions to discrete Riemann-type boundary value problems (RBVP). The convergence 
of s-holomorphic observables is in particular the main tool to establish convergence of Ising interfaces 
to SLE [Smi06, CDHKS12, HoKyll], and to prove conformal invariance of the energy and the spin 
correlations [HoSmlOb, HonlOa, CHI12]. 

A key example is the Ising parafermionic observable of [ChSm09]. On a discrete domain (finite 
simply connected union of faces of Z 2 ), for two midpoints of edges a and z, the observable is defined 



where Z is a partition function, the sum is over collections 7 of dual edges, consisting of loops and 
a path from a to z, and winding (7 : a — > z) is the total turning angle of the path. 

At critical temperature f3 = j3 c , when a is a bottom horizontal boundary edge, the function 
fa '■= f(a, •) is the unique solution of a discrete RBVP: 

• f a is s-holomorphic: for any two incident edges e v = (vu) and e w = (wit), the values of f a 



• On the boundary, values of f a are real multiples of r cw 5 , where r cw is the clockwise tangent 
to the boundary. 

• f a satisfies the normalization condition f a (a) = 1. 

One can then show that the solutions of discrete RBVPs converge to the solutions of continuous 
RBVPs, which are conformally covariant [Smi06, SmilOa, HonlOa, HoSmlOb, ChSm09, Chlzll, 
CHI12]. 

The approach of s-holomorphic functions has proved succesful for the study of conformal invari- 
ance: it applies to arbitrary planar geometries, general graphs and behaves well in the scaling limit. 
Still, the algebraic structures of CFT are not apparent in the s-holomorphic approach: there is 
no Hilbert space of states, no obvious action of the Virasoro algebra and no simple reason for the 
continuous correlations to obey the CFT null-field PDEs. To connect the Ising model with CFT, 
one would like to write algebraic data (e.g. from transfer matrix) in s-holomorphic terms and then 
to pass to the limit. 

1.3. Main results. The goal of this paper is to explore the connection between the transfer matrix 
formalism and s-holomorphicity approach to the critical Ising model, and to lay foundations for a 
quantum field theoretic description that behaves well in the scaling limit. 



by 




satisfy the real-linear equation f a (e v ) + %f a (e v ) = f a (e w ) + \ f a (e w ), where 9 = 



2u—v—w 



\2u—v—w\ ' 
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We construct a discrete analytic continuation matrix, analyze its spectrum and establish a direct 
connection with the Ising transfer matrix. We show that the lattice fermion operators of the trans- 
fer matrix formalism satisfy, as operators, s-holomorphic equations of motion, and we show that 
their correlation functions are s-holomorphic Ising parafermionic observables. Finally, we show that 
Ising correlations can be computed with lattice Poincare-Steklov operators, which encode the geom- 
etry and operator insertions in a manner analogous to the quantum states in the transfer matrix 
formalism. 

The results admit generalizations to non-critical Ising model, with s-holomorphicity replaced by 
a concept of massive s-holomorphicity. 

1.3.1. Discrete analytic continuation and Ising transfer matrix. Let a < b be integers, consider the 
interval I := {x G Z : a < x < b} and let dl := {a, b} denote its boundary. Let I* be the dual of I, 
the set of half-integers between a and b. Write Iq := I* x {0}, Ii := I X {^}, etc. For simplicity of 
notation, we identify edges with their midpoints. 

Lemma (Section 2.4). Let f : Iq — > C be a complex-valued function. Then there is a unique 
s-holomorphic extension h of f to Iq U Ii U 1^ with Riemann boundary values on dli . 

Since s-holomorphicity and RBVP are M-linear concepts, we identify C = R 2 and denote by 
P : CR 2 ) 1 — > (M 2 ) 1 the M-linear linear map / i-> h .In other words, P is the row-to-row 



propagation of s-holomorphic solutions of the Riemann boundary value problem. 

Proposition (Proposition 7 in Section 2.5.). The operator P can be diagonalized and has a positive 
spectrum, given by A^ 1 where \ a > 1 are distinct for a = 1, . . . , |I|*. 

Let / ,: : (C 2 ) 1 — > (C 2 ) 1 be the complexification of P, i.e. the C- linear map such that 
"^ >C I (M 2 ) 1 * = C l^ 2 ) 1 ^ e vector space spanned by the eigenvectors of P c of eigenval- 

ues less than 1 and let P^ : W — > W be the restriction of P c to W . 

Theorem (Theorem 18 in Section 3.3). Let f\ W be the exterior tensor algebra 0^J O A™ W and let 
r (P D C ) : A W -»• A W be defined as 0^' o (P o c ) 0n . Let V + : S + -»• S + be the Ising model transfer 
matrix at the critical point (3 = (3 C , restricted to the subspace S + C S defined as span{e CT : = 1} 
(see Section 1.2.1). 

Then there is an isomorphism p : S + — > A (W ) such that p o V o p^ 1 = Aq x T (Pq ) for some 
A >0. 

It follows in particular that the spectrum of the critical Ising model transfer matrix is completely 
determined by the spectrum of the discrete analytic continuation matrix P. 

1.3.2. Induced rotation and s-holomorphic propagation. The theorem of Section 1.3.1 relies on the 
Kaufman representation of the Ising transfer matrix [Kau49]: V can be constructed from its so-called 
induced rotation Ty on a space of Clifford generators defined below. The connection with discrete 
analysis is made by observing that the s-holomorphic propagation P c is actually equal (up to a 
change of basis) to Ty. 

For k £ I* and a spin configuration a £ {il} 1 - We define the operators p^ : S — > S and 
q k : S ->• S by 

p k (e a ) = cr k+ ie T ( a x for x > k 

i \ ■ 2 , where t x = < 

q k (e a ) = ia k _ie T [-cr x for x < k. 

Let W be the space of operators 5—7-5 spanned by {pk, Qk \ k G I*}. The conjugation O i— >■ VOV^ 1 
defines a linear operator W — > W, which we denote by Ty and call the induced rotation of V. 



DISCRETE HOLOMORPHICITY AND ISING MODEL OPERATOR FORMALISM 



6 



Theorem (Theorem 10 in Section 3.2). Let Ty : W — > W be the induced rotation of V at critical 
point (3 = f3 c , and let P c : (C 2 ) 1 — > (C 2 ) 1 be the complexified s-holomorphic propagation. Then 
there is an isomorphism g : (C 2 ) 1 — > W such that Ty := go P € o g^ 1 . 

1.3.3. Fermion operators. An important tool for the analysis of the Ising model in the transfer matrix 
formalism are the fermion operators; similarly, the study of the scaling limit of the Ising model on 
planar domains relies on s-holomorphic parafermionic observables. We discuss two facts pertaining 
to the relation of these two, namely that the fermion operators are complexified s-holomorphic 
(as matrix-valued functions) and that their correlations are indeed the parafermionic observables 
discussed in Section 1.3.1. 

Theorem (Theorem 19 in Section 4.2). For x £ I*, define the fermion operators ip x -,ipx ■ S S 
by ip x '■= (Px + Qx) and -ip x : = ^ (p x — q x ). Define the operator-valued fermions on horizontal 

edges x + iy G I* x J by ip (x + ty) = V~ y ip x V y and ip (x + iy) = V~ v ip x V y . At the critical point 
f3 = P c , the pair (ip, ip) has a unique operator-valued extension to the edges of I X J, which satisfies 
complexified s-holomorphic equations (see Section If.. 2). 

Conversely, the s-holomorphic parafermionic observables of [SmilOa, ChSm09, HoSmlOb, HonlOa] 
are indeed correlation functions of the fermion operators. 

Theorem (Theorems 22 and 25 in Sections 4.3 and 4.5). The correlation functions of the fermion 
operators are linear combinations of s-holomorphic parafermionic observables. In particular, in the 
box I x J, in the setup of Section 1.2.2, we have 

{iP{z)ip{a)) = f{a,z). 

More generally, all the multi-point correlation functions of ip and ip can be written in terms of 
parafermionic observables. 

This allows one to combine the algebraic content carried by the transfer matrix formalism with 
the analytic content of the s-holomorphicity formalism. As an application we give a simple general 
proof of the Pfaffian formulas for the multi-point parafermionic observables, transparently based on 
the fermionic Wick's formula. 

1.3.4. Operators on Cauchy data spaces. The above results relate the transfer matrix formalism, close 
in spirit to Conformal Field Theory, and s-holomorphicity, suited for scaling limits and conformal 
invariance. We would like to interpret some of the content of the transfer matrix structure in s- 
holomorphic terms. The goal is to pass to the scaling limit and to connect the model with CFT. Can 
we construct quantum states in s-holomorphic terms, that encode domain geometry and insertions, 
and have a scaling limit? 

We present an algebraic construction that encodes the geometry of a domain in a Poincare- 
Steklov operator: all the relevant information about the domain (for correlations) is contained in 
the operator. This operator converges to a bounded singular integral operator in the scaling limit. 

Let Q be a square grid domain with edges £ , let b C d£ be a collection of boundary edges. Let 

_ i 

TZq (resp. Iq) be the Cauchy data space of functions / : b — > C such that / || r CC w on b (resp. 
_ i 

/ || T cw 2 on b), where r ccw = — r cw is the counterclockwise tangent to dQ. 

Lemma (Lemma 28 in Section 5.1). For any u 6 TZq, there exists a unique such that u + v 

_ i 

has an s-holomorphic extension h : — > C satisfying h \\ r cw 2 on 9f2\ b. The mapping u4ti defines 
a real-linear isomorphism Uq : IZq — > Iq. 

The operator Uq is a discrete Riemann Poincare-Steklov operator. The continuous version of this 
operator is defined and studied in [HoPhl2]. 
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When Q = Z x Z + and b = Z x {0}, we have 7£q = IR Z and the operator Uqj, (limit from bounded 
domains) is a discrete analogue of the Hilbert transform (the Hilbert transform maps a function 
w:l->Mtot):l-)-il such that u + v has a holomorphic extension H — > C). 

Proposition (Lemmas 30 and 31 in Section 5.1). The operator Uq is a convolution operator, whose 
convolution kernel is the Ising parafermionic observable at the critical point j3 = j3 c . When 17 = 
I x {0, • • • , N} and b = Io, then Uq is given in terms of the s- holomorphic propagator P N . 

The operators can be used to compute correlation functions by gluing Cauchy data. Denote 
by fa(x,y) the Ising parafermionic observable in domain S7, defined as in Section 1.2.2. 

Theorem (Theorem 35 in Section 5.3). Let fii, 0,2 be two square grid domains with disjoint interiors, 
with edges £±,£2, and let b := d£\ D d£2- The inverse operator Q = (id — U^Uq^ exists. For 
any x G d£\ \ b and any y G £2, the critical Ising parafermionic observable in Q = Qi U SI2 can be 
written as 

fn (x, y) = ^2 foi i k ' x ) Qk,efn 2 {£, y) ■ 
k,eeb 

In other words, the operator Q allows one to 'glue' the domain Q2 to fii, and to compute the 
fermion correlations on Q\ U $^2: all the information about each domain is contained in Uq and 



1.3.5. Away from critical temperature. All the results generalize to temperatures other than the 
critical one. The fermions of Section 1.2.2 satisfy the same boundary conditions and are massive 
s-holomorphic (see Section 2.2 for definition). A massive s-holomorphic propagation Pp : (M 2 ) — > 
(R 2 ) 1 (see Section 2.4) and the non-critical transfer matrix are related like in the critical case. 

Theorem. Let (3 / (3 C . The massive propagator Pp is diagonalizable, with distinct eigenvalues A^ 1 
with X a > 1 for a = 1, 2, . . . , |I*|. 

Theorems of Sections 1.3.1, 1.3.2 and 1.3.3 hold true, if one considers the Ising transfer matrix at 
temperature (3, massive holomorphicity equations, and the massive s-holomorphic propagation matrix 
Pp. 



2. S-HOLOMORPHICITY AND RlEMANN BOUNDARY VALUES 

2.1. S-holomorphicity equations. S-holomorphicity is a notion of discrete holomorphicity for 
complex- valued functions defined on so-called isoradial graphs [ChSmll]. In this paper, we consider 
the case of the square lattice: we consider functions defined on square grid domains, by which we 
mean a finite simply connected union of faces of 1? . More precisely, we will consider functions 
defined on the edges of square grid domains; when necessary, we will identify these edges with their 
midpoints. 

S-holomorphicity is a real-linear condition on the values of a function at incident edges; it implies 
classical discrete holomorphicity (i.e. lattice Cauchy-Riemann equations) but is strictly stronger. 
The fact that Ising model parafermionic observables are s-holomorphic is the key to establish their 
convergence in the scaling limit and hence to prove conformal invariance results. 

Definition 1. Let $7 be a square grid domain. Set A := e l7r / 4 . We say that F : $7 — > C is s- 
holomorphic if for any face of £1 with edges E, N, W, S (see Figure 2.1), the following s-holomorphicity 
equations hold 
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FIGURE 2.1. The four edges adjacent to a face on the square lattice. 



(2.1) F(N) + XF(N) = F(E) + XF(E) 

F(N) + A _1 F(iV) = F(W) + X~ l F(W) 
F(S) + A 3 F(S) = F(E) + X 3 F~(E) 
F{S) + A~ 3 F(S) = + A' 3 F(W0. 

In other words F: £1 — > C is s-holomorphic if for any pair of incident edges e v = (uv) and 
e w = (uw), we have F (e v ) + %F (e v ) = F (e w ) + ^F(e w ), where 9 = ^Z V V ZZ\ ■ Equivalently, the 

orthogonal projections (in the complex plane) of F (e v ) and F (e w ) on the line coincide. 

The above equations imply (but are not equivalent to) the usual lattice Cauchy-Riemann equa- 
tions: for the four edges around a face as in Figure 2.1 we have F (TV) — F (S) = \ (F (E) — F (W)), 
and a similar equation holds for the four edges incident to a vertex. Discrete Cauchy-Riemann 
equations imply in turn the discrete Laplace equation Ylz=x±i x±i — F (X)) = for every 

edge x en\ on. 



2.2. Massive s-holomorphicity. We now define a perturbation of s-holomorphicity which we call 
massive s-holomorphicity. The massive s-holomorphicity equations with parameter (3 are M-linear 
equations satisfied by the Ising model parafermionic observables at inverse temperature (5. At the 
critical point f3 = f3 c , massive s-holomorphicity reduces to s-holomorphicity. 

Definition 2. Let f3 > 0, let v = v (/3) be the unit complex number be defined by v = A ^z^, where 

a = e~ 2 @ and A = e i7r//4 . A function F : £1 — > C is said to be massive s-holomorphic with parameter 
(3 if for any face of $7 with edges E, N, W, S, we have 



(2.2) 



F(N) + v^XF^N) 
F(N) + vX^FjN) 



-i 



F(E) + XF(E) 



F(S) + vX 3 F(S) 
F(S) + u~ 1 X~ 3 FjS) 



vF{W) + X^FiW) 
vF(E) + X 3 F(E) 
v~ 1 F{W) + X~ 3 F(W). 



At /3 = f3 c , we have v = 1 and these equations coincide with the Equations (2.1) defining s- 
holomorphicity. It can be shown (see [BeDCIO]) that massive s-holomorphicity implies (but is not 
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equivalent to) the massive Laplace equation 

\ £ (f(z)-f(x)) = »f(x) vxen\on, 

Z=X±l,X±i 

with the mass fi = fi((3) given by fx = — 1, where S = sinh (2j3). The dual inverse temperatures 

/? and /3*, related by sinh (2/3) sinh (2/3*) = 1, have equal masses //(/?) = fi(f3*), and at the critical 
point the mass vanishes Li((3 c ) = 0. 

2.3. Riemann boundary values. The boundary conditions that are relevant for the study of Ising 
model specify the argument of a function on the boundary edges d£l: these conditions are trivially 
satisfied by the Ising parafermionic observables for topological reasons (see Section 4.3), at any 
temperature. 

Let Q be a discrete square grid domain. The boundary d£l of f2 is a simple closed curve. For an 
edge e G d£l, this defines a clockwise orientation r cw (e) of e, which we view as a complex number: 
t cw (e) G {±1} if e is horizontal and r cw (e) = {±i} if e is vertical. 

_ i 

Definition 3. We say that a function / : — > C satisfies Riemann boundary conditions / || r cw 2 
at an edge z G dQ if 

_ i 

f(z) || r cw 2 (z) , 

_ i 

i.e. / (z) is a real multiple of r cw 2 (^)- 

_ i 

When Q is a rectangular box I x J, the condition / || r cw 2 means that / is purely real on the top 
side of Q, purely imaginary on the bottom side, a real multiple of A = e l7r//4 on the left side and a 
real multiple of A -1 = e~ l7r / 4 on the right side. 

2.4. S-holomorphic continuation operator. For a (massive) s-holomorphic function on a rect- 

_i 

angular box with Riemann boundary conditions || r cw 2 , we can propagate its values row by row as 
illustrated in Figure 2.2. This is supplied by the following lemma (we use the same notation as in 
Section 1.3.1). 

Lemma 4. Consider the box Ix {0, 1} for an integer interval I = [a, b]nZ and let I* = [a, 6]n(Z+ \) 
be its dual. 

Let f : Iq — > C be a complex-valued function and let f3 > 0. Then there is a unique massive 
s-holomorphic extension h of f to Iq U Ii U I* with Riemann boundary values on dli . 

Proof. For z G Ii \ dl i , the value h (z) can be solved uniquely from the last two of Equations (2.2) 

in terms of / (z — \ — §) and / {z + \ — §)■ For z G dl i , the value h (z) can be solved uniquely 

from the Riemann boundary condition and (2.2) in terms the value / (2 — | ± 5) (± depending on 
whether z is on the left or the right part of dli). For z €.!]_, h(z) can be solved in terms of h at 

z + 7} — 5 Gli and 2—3 — 3 e ^- tne ^ rst two °^ Equations (2.2). The definition of h thus 
obtained satisfies all the required equations. □ 

Definition 5. Let I be an interval of Z as above. We define the /3-massive s-holomorphic propagator 
Pp : (M 2 ) 1 —t- (M 2 ) 1 by Ppf = where / and h are as in Lemma 4 

We can explicitly write down the s-holomorphic propagator in critical and massive cases. The 
explicit form will be useful in the next section. 
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FIGURE 2.2. The values of a massive s-holomorphic function in the row above can 
be solved in terms of the values in the row below, both in the bulk and on the 
boundary. 



Lemma 6. Let I* = {a + 5, a + §,...,& — ^} and denote the left and right extremities by ki = a 
and k R = b — \. Set A := e l7r//4 . The s-holomorphic propagator P is given by 



' (P f) (k) = f(k-l) + 2 + $ /( fc + 1) 

+ 71 ^ k " !) - /(*0 + ^ /(* + 1) Vfc G I* \ {k L , k R } 

(Pf)(k L )= (l + ^f(k L ) + ^f(k L + l) 

+ (a 3 + ^)7M+^7(^TT) 

(Pf)(k R )= ^f(k R -l)+(l + ±)f(k R ) 

+^7(^T) + (a- 3 + ^)7(M. 



For /3 7^ /3 C; denote S := sinh(2/3) ; C := cosh (2/?). TTie massive s-holomorphic propagator Pp 
given by 
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' (P fi f) (k) = =§j *f(k-l) + ff(_k) + =^f(k + 1) 

+gf{k - 1) - cf(k) + gf{k + i) vfc € r \ {k L , k R } 

{Ppf)(k L )= £g£ /(* L ) + =g* /(fc £ + 1) 

i - (y y- si /faj+i/feTij 
(p/, /)(*«)= =^ - 1) + (^f) 

+S7(i^+( - (fl nr c - s) )7M- 

2.5. Spectral splitting of the propagator. 

Proposition 7. T/ie matrix P 13 is symmetric, with eigenvalues A^ 1 , where X a > 1 are distinct for 
a = 1, . . . , |I|*. 

Proof. Clearly P^ is invertible: the inverse of Pg is the propagation of values of a massive s- 
holomorphic function downwards. Notice that exchanging S and N in the massive s-holomorphic 
equations (2.2) amounts to replacing F by iF. Denoting by j : (M 2 ) 1 — > (M 2 ) 1 the (real-linear) 
involution / i-> i/, we deduce that (P/g) -1 = j ° Pp ° J -1 . We deduce that the spectrum of Pg is 
the same as the one of (Pp)~ and hence that the eigenvalues are of the form A^ 1 with \ a / for 
ol = 1 , . . . , 1 1* | . 

For i) £ C, observe that the real-linear transpose of the map z ^ i|z is z H> p and that the 
map z i — y rfz is real-symmetric. From the formula of Lemma 6, we deduce that Pp is symmetric. 
To show that Pp is positive definite, it is enough to show this at the critical temperature f3 = (3 C , 
since the eigenvalues of P^g are continuous in (3 and cannot be zero. We can write the propagator 

Pp as BA, where A : (R 2 ) 1 — > (M 2 ) 1 is the propagation of / : Iq — > C to Ii (see Definition 2) and 
B : (R 2 ) 1 -> (M 2 ) 1 * is the propagation of g : Ii ->• C to T[. At (3 = (3 C , we have that B = A T and 

hence P = A T A is positive definite. 

Let us now show that 1 cannot be an eigenvalue. Suppose / : Iq — > C is such that Ppf = /; we 
want to show that / = 0. Let h : Iq U Ii U 1^ — >■ C denote the massive s-holomorphic extension of /. 

At (3 = (3 C (i.e. when Pg = P), there is a particularly simple argument. Since h satisfies the 
discrete Cauchy-Riemann equations, i.e. for any x G I * , we have 

2 

h { x + 2) " h { x ~ 2) = \{ h { x + 2) " h ( x " 2) ) = T (/ {x) " J {x)) = a 

Hence /i must be constant on Ii and the Riemann boundary conditions easily imply that h = 0. In 

2 

turn, this implies that / = 0, by the s-holomorphicity equations. 

For general (3, writing Ii = {xl,xl + 1, . . . , Xr}, we can deduce (from an explicit computation) 
that h (x + 1) and h (x) (for x £ {xl, ■ ■ ■ , xr — 1}) satisfy a linear relation: 

(2.3) h(x + l) + \Bh {x + 1) = h(x)- \Bh (x), 

for S = v/^Qmfi/), where z/ = e i3,r / 4 ^±i anc i a = e ~ 2 P as above. 

The Riemann boundary condition on the left extremity imposes that h(x L ) = Ce" i7r / 4 for some 
Cel, and the above equation (2.3) yields h(x) = Ce~ i7T / 4 A x ~ XL where A = j^f 7^ 0. But the 
Riemann boundary condition on the right extremity, /i (x/?) || e i7r / 4 , then requires that CA XR ~ XL = 
and hence h = everywhere by the massive s-holomorphic equations. 

Finally we show that the eigenvalues are distinct. Suppose A > is an eigenvalue of Pp, and 
let f\ £ (M 2 ) 1 * is an eigenvector, and let h\ be the massive s-holomorphic extension of f\ to 
Ig U Ii U I*. The massive s-holomorphicity equations can be solved to obtain a recursion relation 
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h\(x + 1) = r}h\{x) + rf'h\(x) with some explicit 77,77' £ C. This shows that the eigenspace is 
one-dimensional. □ 

3. Transfer Matrix, Clifford Algebra and Induced Rotation 

In this section we review fundamental algebraic structures underlying the transfer matrix formal- 
ism introduced in [Kau49]. See [Pal07] for a recent exposition with more details. 

3.1. Transfer matrix and Clifford algebra. In the introduction, Section 1.2.1, we defined the 
Ising transfer matrix V : S — > S with fixed boundary conditions at the two extremities of the row 
I = {a, a + 1, . . . , b — 1, 6} as the product 



V = (V h )z V v (V h )^ 

where the matrix elements in the basis (e CT ) indexed by spin configurations in a row, a G {il} 1 , are 
given by 

r exp (j3 Y? i=a <?iP?J if cr a = p a and a h = p b 
otherwise 



and 



^ Up ' [0 otherwise. 



There is a two-fold degeneracy in the spectrum of the transfer matrix: the global spin flip a H> —a 

commutes with both V v and (^ h ) 5 . To disregard the corresponding multiplicity of eigenvalues, we 
restrict our attention to the subspace 

S + = span {e a | cr fe = +l} , 

spanned by spin configurations that have a plus spin on the right extremity of the row. This subspace 

1 

is invariant for both V v and (V h ) 2 . Note that the dimension is given by 

dim(cS+) = 2l J l- 1 = 2l p l = 2 b ~ a . 
In Section 1.3.2, we defined also the operators pk and % on 5, for A; G I* = {a+^, a + |, b — |, b — j}, 



by 



Pk (e CT ) = ^ + ie T j a x for x > k 

1 \ ■ , where t x = < 

q k (e a ) = io- fc _ie r for x < k. 



Again, the subspace S + C S is invariant for all Pk,qk- We denote W = spanjp^.,^ | k G I*} C 
End(5). 

It is easy to check that p^, qt satisfy the relations 

PkPt + PiPk = 2 5 k>e ids 
Qkqe + qeqk = 2 S k ,e id s 
PkQe + qiPk = o, 

i.e. that they form a Clifford algebra representation on <S_|_ and on S. This representation is faithful, 
so we think of the Clifford algebra (£tiff simply as the algebra of linear operators S — > S generated 
byW. 
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Consider the symmetric bilinear form (-, •) on W given by (pk,Pi) = ^k,u (<?fc> Ql) = ^k,h (Pk, Ql) = 
0. Then Cliff is the algebra with set of generators W and relations uv + vu = (u, v) 1, for u, v £ W. 
The dimensions of the set of Clifford generators and the Clifford algebra are 

dim(W) = 2|I*| = 2(6 - a), dim(£liff) = 2 dim ( w ) = 2 2 I F I = 2 2 ^ a ). 

The transfer matrix can be written in terms of exponentials of quadratic expressions in the Clifford 
algebra generators as follows. 



Proposition 8. We have 



(Vf = exp (if E 



QkPk 



feel* 



V v = e^(2sfl^ expli/3* £ p.^q. 

\ jei\di 

where (3* is the dual inverse temperature given by tanh (/?*) = e~ 2/3 and S = sinh (2/3). 
Proof. Note that the operator ig/j^ has the following diagonal action in the basis (e CT ) 

iqkPk e CT = a k+ ia k _i e a , 

so the first asserted result e l ^^- /kQkPk = (V^ 1 )? follows immediately. The operator ip,-_iO,-, i inverts 
the value of the spin at j, 



-5^ + 5 



tP 7 -_i9,-+i e CT = e a >, where ^ = < 

J I — cjj for x = j, 

so we have 

exp (i/3V_i<? 1+ i ) e a = cosh(/3*) e CT + sinh(/3*) e a >. 



a x for x / j 



Taking (3* such that tanh (/?*) = e 2/3 and computing the product of these commuting operators at 
different j we get 

ji-i \ 

3(J = cosM/?*) 6 -"- 1 tanh(/3*) # ^ CT ^"^e p 




pefii} 1 

Pa=&a, Pb = <J b 



e -2/3 ' 



6-1 



(25) 1 -V £ exp £ 

CiPi e p , 

pefil} 1 \i=a+l / 

Pa=0"a, Pb = (J b 

which is the second asserted result. □ 

3.2. Induced rotation. Since the constituents of the transfer matrix are exponentials of second 
order polynomials in the Clifford generators, conjugation by the transfer matrix stabilizes the set of 
Clifford generators. In the formulas below, we use the notation 

s = sinh(/3) s* = sinh(/T) 

c = cosh(/3) c* = cosh(/T) 

S = sinh(2/3) = . * - C = cosh(2/3) - cosh (W 



sinh(2/3*) " v ' ' sinh(2/3*) ' 

The following lemma is a result of straightforward calculations, which can be found e.g. in [Pal07]. 
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Lemma. Conjugation by V h is given by the following formulas on Clifford generators Pk,qk (k £ 1*) 



i ,i 



V h J 2 op k o (V h y = cp k -isq k 

(v h y*oq k o (v h Y =isp k + cq k . 

Let ki = a + \ and k R = b — | be the leftmost and rightmost points of I*. Conjugation by V v is 
given by 

C i 

(y v ) _1 o Pk oV v = —p k + - q k+ i for k^k R 
— i C 

(V^ 1 o q k o V Y = — + — q k for k / k L 

and on the remaining generators by 

(VT 1 o PkR o V v = p kR (O -1 Qk L o y v = % L • 

We see that W C End(5) is an invariant subspace for the conjugation by the transfer matrix V. 
The conjugation is called the induced rotation of V, and denoted by 

Ty.W^W T v (w) = VowoV' 1 . 

Note that the induced rotation Ty preserves the bilinear form, (Tyu, Tyv) = (u, v) for all u, v G W. 

We will next show that the induced rotation is, up to a change of basis, the complexification of 
the row-to-row propagation Pp of massive s-holomorphic functions satisfying the Riemann bound- 
ary condition. To facilitate the calculations, we introduce two symmetry operations on the set of 
Clifford algebra generators. We define a (complex) linear isomorphism R : W — > W and a (complex) 
conjugate-linear isomorphism J : W — > W by the formulas 

R{Pk) = iq a +b-k R{qk) = ~ iPa+b-k (extended linearly) 

J(pk)=ipk J(qk) = —iqk (extended conjugate-linearly). 

We will moreover use for W the basis 

i>k = -^(pk + qk), 4>k = ^(pk- qk), fee I*. 

Lemma 9. The maps R and J commute with Ty, i.e. we have 

Ty o R = Ro Ty, Ty O J = J O Ty. 

For all k G I* we have 

R(lpk) = 4>b+a-k, R(tpk) = lpb+a-k 

J{i>k) = i>k, J{4>k) = iik- 

Proof. By the explicit expressions of Lemma 3.2 one easily verifies that R and J commute with the 

I 

conjugation by both (V h ) 2 and V v . □ 

Theorem 10. The induced rotation Ty : W — > W is up to a change of basis equal to the complexi- 
fication of Pp, i.e. there exists a linear isomorphism g : (C 2 ) 1 — > W such that Ty = g o P c o g -1 . 

Proof. Consider the action of the induced rotation in the basis (V , fc>V'fc)fc€l*- With the formulas of 
Lemma 3.2, it is straighforward to compute that for feel*\ (<9I*) 

T v \^k) = C s^+{-\- ^ ^-i + ("2 + 2^) ^ 

nl 1C 7 lC 7 

- Ci^k + 2^^-i + 2^Wi' 
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To get a formula for Ty 1 ^^), k G I*\(<91*), apply the map J on this and use Lemma 9. We still need 
formulas for the two extremities, k € dl* = {a + ^, b — |}. On the left extremity, at ki, = a + a 
straightforward calculation yields 

rr-i,, , C(S + C) (1 + xS). 



-5(C + S)+i(C-S) 7 _c 
25 ^ + 2S 



To get a formula for Ty (i(>k L ), apply the map J on this. To get the formula for Ty (Y^), where 
kn = b — \, apply the map R. To get a formula for Ty 1 {ipk R ), a Pphy the composition J o R. 

Since the coefficients in the formulas for Ty 1 (ipk) coincide with the coefficients in the formulas for 
(Ppf)(k) in Section 2.4, and the coefficients in the formulas for Ty 1 ('4>k) are the complex conjugates 
of the corresponding ones, we get that the complexification of Pp agrees with Ty 1 up to a change of 
basis. This finishes the proof, because Pp and its inverse P^ 1 are conjugates by Proposition 7. □ 

3.3. Fock representations. Finite dimensional irreducible representations of the Clifford algebra 
Cliff are Fock representations, defined below. To define a Fock representation, one first chooses a 
way to split the set W of Clifford algebra generators to creation and annihilation operators. Let 
(•, •) denote the bilinear form on W defined in Section 3.1. A polarization (an isotropic splitting) is 
a choice of two complementary subspaces Wm (creation operators) and Wb (annihilation operators) 
of the set of Clifford algebra generators, 

W = W ffl 9W B , 

such that 

(wa, w' m ) = for all u> ffl , w' m £ Wm 
(wb, w' b ) = for all u>b, w' b € Wb. 

Note that due to the nondegeneracy of the bilinear form (•,•), the two subspaces Wb and Wb are 
naturally dual to each other, and in particular 

dim(W ffl ) = dim(W B ) = ^dim(W) = |I* | = b - a 

is the number of linearly independent creation operators. 

As a vector space, the Fock representation corresponding to the polarization W = Wb © Wb is 
the exterior algebra of Wb , 

|i*| 

A w ffl= 0( A "Wffl)- 

n=0 

To define the representation of the Clifford algebra on this vector space, let (fia)aJi he a basis of 
Wb and (oq)^! the dual basis of Wb, i.e. {o) a ,ap) = 5 a ^. The action of the Clifford algebra on 
the Fock space f\ Wb is given by the linear extension of the formulas 

A a h A • • • A a L) = a « A a k A a k A • • • A a L 

n 

a a k a ■ ■ ■ a a k) = j2^- i y~ l6 ^ a k a • • • a a k-i a a • • • a °l- 
3=1 

The vector 1 G C = A°W ffl C A W ffl is called the vacuum of the Fock representation: it is annihilated 
by all of Wb- Irreducible representations are characterized by such vacuum vectors as follows. 
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Lemma 11. Suppose that W = Wm © Wb is a polarization. Any irreducible representation of Cliff 
is isomorphic to the Fock representation f\Wm- If a representation V of Cliff contains a non-zero 
vector t> vac G V satisfying Wa^vac = 0, then the Fock space /\ Wm embeds in V by the mapping 

°k A a k A • • • A °L ^ (°k a k " " " 4J- V ™- 

Proof. Consider a representation V of Cliff. Choose a non-zero v ^ G V. Define recursively t/ a ), for 
a = 1,2, . . . |I*|, to be a a v^ a ^^ if this is non-vanishing and tX" -1 ) otherwise. Then f vac = t/I 1 *!) is 
non-zero and We^vac = 0, i.e. v vac is a vacuum vector. This argument shows in particular that any 
non-zero subrepresentation of the Fock representation /\ contains the vacuum vector 1 G A°Wh, 
and thus the Fock representation is irreducible. The mapping f\ Wm — > V given in the statement 
defines a non-zero intertwining map of Clifford algebra representations, and by irreducibility of the 
Fock representation, this is an embedding. If V is irreducible the embedding must be surjective, and 
thus an isomorphism. □ 

The fact that the Fock representation is the only isomorphism type of irreducible representations 
of the Clifford algebra would follow already from the irreducibility of the Fock representation and 
the observation that (dim(/\ Wa)) 2 = dim(Cliff), by the (Artin-)Wedderburn structure theorem. 

A standard tool for performing calculations in the Fock representation is the following. We recall 
that the Pfaffian Pf(-A) of an antisymmetric matrix A G C nxn is zero if n is odd, and if n = 2m is 
even, then it is given by 



1 m 



^ M2o-l),ir(2 a )- 
Tre&m s =l 



Lemma 12 (Fermionic Wick's formula). Let W = We]©VVb be a polarization, and consider the Fock 
representation f\Wm- Letv vac = 1 G A°Ws] C /\Wa be the vacuum andv* ac G (a°Ws])* C 
be the dual vacuum normalized by (v * ac , v vac ) = 1. Then for any <f>i, . . . , (j) n G W we have 



«ac> 01 ■ ■ ■ ^n^vac) = Pf ^ « ac , <f>i4>jV vac ) 



Proof. Write the elements <fii G W as sums of creation and annihilation operators, and then anti- 
commute the annihilation operators to the right and the creation operators to the left. □ 

3.4. A simple polarization for low temperature expansions. The following lemma gives one 
of the simplest possible polarizations. 

Lemma 13. The following formulas define a polarization 

w4 +) = span [p k - \q k \ k G I* } 

Vti +) = span {p k + xq k \ k G I*} . 

Proof. The vectors p k + iqk, Pk — Wk f° r k G I* form a basis of W, and we have (p k ± iq k , pi ±iqi) = 
25 M + + 0+(±i) 2 2,5 M = 0. ' ' ' □ 

Recall that in the state space S of the transfer matrix formalism we have the vectors corresponding 
to the constant spin configurations in a row, 

e ( +)GcS+ (+) = (+1,+1,...,+1) G {il} 1 

e ( _)G<S (-) = (-l,-l,...,-l)€{±l} 1 . 

Directly from the defining formulas of the operators Pk,Qki one sees that the vectors e( + ),e(_) G S 
satisfy (p k + \q k )e^ = and (p k + i<Zfc)e(_) = for all k G I*. 
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Corollary 14. As a representation of the Clifford algebra, 5+ is isomorphic to the Fock space 

Aw^, with vacuum vector vitc = e (+), and S is isomorphic to the direct sum of two copies of this 
Fock space. 

We emphasize that the polarization of this subsection is not the physical one, but by Lemma 
11, the isomorphism type of the Fock representation doesn't depend on the polarization, so the 
state space of the transfer matrix formalism is in fact a Fock representation for any polarization. 
The polarization is, however, the zero temperature limit (/3 /* oo) of the physical polarizations 
of the next section, and it is very closely related to the low temperature graphical expansions of 
correlation functions and observables considered in Sections 4.3, 4.4 and 4.5. In particular, a slight 
modification of this simple polarization together with the fermionic Wick's formula will be used for 
the proof of Pfaffian formulas for fermion operator multi-point correlation functions and multi-point 
parafermionic observables. The modified polarization is the following. 

Lemma 15. Let N € N. The following formulas define a polarization 

W^ y ' N = sp a n{V- N (p k - iq k )V N \ k € I* } 

= span {p k + iq k \ k G I*} 

for all f3 except possibly for isolated values. The space 5+ is isomorphic to a Fock representa- 



tion /\ (w£ +);JV ), w % th vacuum vector vitc' N = e (+) £ <5+ and dual vacuum vector (v 



vac j 



e 



(+) v e (+) 

Proof. The special case N = was treated above. Since we have {p k ± \q k , p\ ± iqi) = and the 
bilinear form is invariant under Ty, (Ty N u, Ty N v) = (u,v), it follows that also for general N the 
choice of subspaces is a polarization if the two subspaces span the whole space W, that is if the 
vectors Ty N (p k — \q k ) and p k + \q k for k G I* form a basis of W. It suffices to show that the 

matrix (p k + iq k , Ty N (pi — iqi)) of the bilinear form is non-degenerate. The non-degeneracy 

J fc,zei* 

is evident in the limit (3 /* oo, since e~ 2 @Ty 1 (p k — iq k ) = Pk~ ^Qk + 0(e~P) and e~ 2 ^Ty 1 (p k + iq k ) = 

0{e~P) by the formulas of Lemma 3.2. The determinant det I (p k +iq k , Ty N (pi — iqi)) 

is analytic in (3, so its (possible) zeroes can't have accumulation points. We conclude that W = 
is a polarization except possibly for isolated values of (3. 



The same calcuation as before shows that e( + ) is a vacuum vector of the Fock representation 

"+)' 

(+)' 



S+ = f\ W^' N , and similarly from the calculation ej + ^(p k — iq k ) = we get that the dual vacuum 
of S+ A W^ )]N is proportional to eJ.,V N . □ 



3.5. The physical polarization. The relevant polarization and basis is the one in which the 
particle-states oj^aj^ • • -a^.f vac are eigenvectors of the evolution defined by the transfer matrix. 
We make use of the fact that 1 is not an eigenvalue of Ty, which follows from Proposition 7 and 
Theorem 10. 

Lemma 16. Let W^ hys C W be the subspace spanned by eigenvectors of Ty with eigenvalues less 
than one and Wg hys C W the subspace spanned by eigenvectors of Ty with eigenvalues greater than 
one. ThenW^ ys @W^ ys is a polarization. As a representation of Cliff, the space 5+ is isomorphic 
to the Fock representation f\yV^ ys . 

Proof. Recall that for any u, v £ W we have (Tyu, Tyv) = (u,v). For eigenvectors of u,v of Ty it 
follows that (u, v) can be non-zero only if the eigenvalues are inverses of each other, and thus the 
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bilinear form vanishes when restricted to Wg hys or yV^ ys . Finally, W = W^ hys ® Wg hys because Ty 
is diagonalizable with real eigenvalues and 1 is not an eigenvalue. □ 

Then let (Ottjfjj be a basis of Wg hys consisting of eigenvectors of the induced rotation Ty{a a ) = 
X a a a with X a > 1, and let (aa)|JJ 1 be the dual basis of W^ hys , i.e. (aq,^) = 5 a ^. Note that we 
have Ty(ai) = A^aL 

Proposition 17. If v £ S is an eigenvector of V with eigenvalue A, then the vector alt) G S is 
either zero or an eigenvector with eigenvalue A Q J 1 A and a a v £ S is either zero or an eigenvector of 
eigenvalue X a A. In particular, if Aq is the largest eigenvalue ofV and t>va yS G 5+ is the corresponding 
eigenvector, then t>va ys is a vacuum of the Fock space S+ and the vectors o} ai o) a2 ■ ■ ■ a>a n .v vac form a 
basis of S + consisting of eigenvectors with eigenvalues Aq x n"=i 1 ■ 

Proof. For v G S an eigenvector, Vv = Av, compute Va a v = (VaaV^^Vv = Ty(a a )Vv = X a Aa a v, 
and similarly for a) a v. It is then clear that v^aF is annihilated by all of Wb, because A a Ao is larger 
than the largest eigenvalue of V. □ 

Theorem 18. LetP^: (C 2 ) 1 * ->• (C 2 ) 1 * be the complexified massive s-holomorphic row-to-row prop- 
agation, and let W C (C 2 ) 1 be the subspace spanned by eigenvectors of Pp with eigenvalues less 

than one. On the exterior algebra f\W Q = 0^' o A n W define T(Pj~) = ©^(^IvKo) ™- Then 
there is a linear isomorphism p : S+ — > f\ W such that 

p o V o p- 1 = const, x r(jf). 

Proof. The state space 5+ is isomorphic to the Fock space f\ Wg hy& by Lemma 16. By Proposition 

17, in this identification the transfer matrix V becomes diagonal in the basis A • • • A da n , with 

11*1 / \®n 

eigenvalues Ao Ils=i S an< ^ tnus ^ coincides with 0„ =o ( Tv\ w ph ye j apart from the overall 
multiplicative constant Ao It remains to note that by Theorem 10 the induced rotation Ty. W — > 
W coincides up to isomorphism with the complexification P^ : (C 2 ) 1 * — > (C 2 ) 1 * of the row-to- 
row propagation, and the same holds for the restrictions Tj^phys and Pp\w a to the corresponding 
subspaces. □ 



4. Operator Correlations and Observables 

In this section we discuss correlation functions of operators in the transfer matrix formalism. 
We introduce in particular holomorphic and antiholomorphic fermion operators, and show that they 
form an operator valued complexified s-holomorphic function. The low temperature expansions of the 
fermion operator correlation functions are simply expressible in terms of parafermionic observables. 



4.1. Operator insertions in the Ising model transfer matrix formalism. We consider the 
Ising model in the rectangle I x J, with I = {a, a + 1, . . . , b — 1, 6} and J = {0, 1, . . . , N — 1, N}, and 
we denote the row y by I y = I x {y}. We use the notation of Section 3.1 for the transfer matrix (with 
locally constant boundary conditions on the left and right sides of the rectangle) and the Clifford 
algebra. 

The total energy of a spin configuration s G {±l} IxJ is H(s) = -Ed~)c s « s «" with the sum 
over v,w G I x J that are nearest neighbors on the square lattice, \v — w\ = 1. The probability 
measure of the Ising model with plus boundary conditions is given by P + [{s}] = -^+e~P H< ^ on 

the set js G {±l} IxJ s|a( IxJ ) = +1 j of spin configurations that are +1 on the boundary of the 
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rectangle. The normalizing constant in the formula is the partition function 

z+ = y, e ~ ms) - 

se{±i} IxJ 
s b(ixJ)=+ 1 

The partition function can be expressed in terms of the transfer matrix V by expanding a product 
of transfer matrices in the basis (e CT ) indexed by spin configurations in a row, a £ {±1}\ More 
precisely, we have 

Z+ = Yl e^^ s ^ = i T V N i =: (i\V N \i), 
se{±i} IxJ 

s ls(IxJ) = + ! 

where the "initial state" i and the "final state" f are given by i = f = (V h ) 5 e( + ) = e 2 ^ 'e( + ) — we 
included a factor to correctly take into account the interactions along the horizontal edges in the 
top and bottom rows. 

The spin operators dj : S — > S are the diagonal matrices in the basis (e CT ) with diagonal entries 
given by the value of a G {il} 1 at position j G I, i.e. 

&j(e a ) = Uj e a . 

Note that for example the expected value of the spin s z at z = x + iy G I x J, with respect to the 
probability measure P + of the Ising model with plus boundary conditions, can be written as 

F+ r , = E s ^exp(/?E^^s w ) = (f\v N -ya x yy\i) 

[ZI E s exp(/3E^s,s w ) <f|^|i> ' 

by expanding also matrix products in the numerator in the basis (e CT ). Moreover, the initial and 
final states i, f oc e( + ) could be replaced by e( + ) because the constants would cancel in the ratio. 
Finally, the formula takes a yet simpler form if we define the time-dependent spin operator 

&{x + iy) = V~ y a x V y 

and indeed it is simple to check that then 

_ (e(+)|^ JV 'r(2i)---'T(2r)|e( + )) 



.1=1 



(e (+) |^|e (+) ) 



We define, as in the proof of Theorem 10, the Clifford algebra elements ipk = ~^(Pk + Qk) £ VV, 
4>k = ^(Pk — Qk) S W for S I*. The corresponding time-dependent operators 

il>(k + iy) = V-y^ k Vy ' K€L,y€J, 

are called the holomorphic fermion and the anti-holomorphic fermion, respectively. The reason 
for this terminology is Theorem 19 below, which states that the pair of operator valued functions 
(if), ■0) satisfies local linear relations that have the same coefficients as the defining relations of s- 
holomorphicity for a function and its complex conjugate. 
The following abbreviated notation 

_ <e (+) |VV%i)---V> (n ^n)|e (+) > 



(V (1) (*l)"-^ (n) (*n))J 



'ixJ • (e (+) |^|e (+) ) 

will be used for the correlation functions of the fermion operators, where z\, . . . , z n are edges, and 
for each i = 1, 2, ... n we let ip^ stand for either ip or ip. 

Note that if R: W — > W is the linear isomorphism introduced in Section 3.2, then we have 

R^(z))=i,(r(z)), R^(z)) = ^(r(z)), 
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where r(x + \y) = a + b — x + \y, and if J : W — > W is the conjugate-linear isomorphism of the same 
section, then 

jy>(z)) = $(z), j$(z)) = i>(z). 

4.2. S-holomorphicity of fermion operator. The fermion operators tp(z) and tp(z) were defined 
in the previous section for z G I* x J, i.e. on the set of horizontal edges of the rectangle I x J. The 
following theorem says that we can extend to vertical edges so that the pair (tp, tp) is a complexified 
operator valued (massive) s-holomorphic function. 

Theorem 19. Let the fermion operators tp(z),tp(z) be defined by Equation (4-1) for horizontal edges 
z G I* x J. Then there exists a unique extension of tp and tp to the set of vertical edges I x J* , such 
that the following local relations hold. For any face, with E, N, W, S the four edges around the face 
as in Figure 2.1, we have 



(4.2) i/,(N) + v^X^N) = i/- V(£) + \$(E) 

il>(N) + v\- l 4>{N) = wp(W) + A _1 ^(W) 

il>(S) + uX 3 tf;(S) = utP{E) + X 3 $(E) 

tp(s) + u^x'^is) = i/"V(w) + x~ 3 ij{w), 

and on the left and right boundaries we have 

(4.3) tp(a + W) + i i>{a + \y) = 

ip(b + \y) — \ V>(6 + iy) = for any y G J* . 

Remark. The coefficients of the linear relations among the operators on incident edges, Equations 
(4.2), coincide with the coefficients in the definition of massive s-holomorphicity, Definition 2. Simi- 
larly, coefficients in the Equations (4.3) coincide with the equations defining the Riemann boundary 
condition on the left and right boundaries. The situation at the top and bottom boundaries is 
slightly different: the operators -ip and -p are linearly independent, but when the operators are ap- 
plied to specific boundary states we recover similar relations, e.g. at the bottom for z G Io we have 
(^(z) + #s))e (+) =0. 

Proof. The uniqueness of such extension is clear by the following explicit construction similar to 
the one in the proof of Lemma 4. Consider the vertical position y G J*. For z G (1^ \ dl y ) one 
can solve for ip(z) from Equations (4.2) (the third and fourth equations on the plaquettes on the 

left and right of z) in terms of the operators ifj(w) and 4>(w), w G I*_i, more precisely in terms of 

y 2 

ip(z— \ — 5), ip{ z ~ \ ~ 5)) 4 , { z ~\~ \ ~ 2)' 4 , { z ~\~ \ ~ 3)- Similarly one can solve for ip(z) and the result 
is J(ip(z)). For z G dl y = {a + iy, b + \y} on the boundary, using both Equations (4.2) and (4.3) 
one can solve for ip{z) in terms of the operators ip(w) and tp(w), w £l* l , more precisely in terms 

of ip(a+ 5 — 5) and ip(a+ \ — ^) or ip{b— 3 — 5) and ip(b— \ — ^). Again similarly ip(z) = J(ip(z)). 
Extending ip and tp to I y with the above formulas in terms of tp and tp in the row I*_i, the 

y 2 

Equations (4.3) as well as the third and fourth of Equations (4.2) hold by definition. Then note that 
by a similar argument, there are unique values of tp and tp in the row I* + i such that the first and 

second of Equations (4.2) hold. Since the coefficients of the equations we have used are the same 
as the coefficients defining massive s-holomorphicity, the unique definitions of tp in the row I* 1 are 

expressible as linear combinations of the operators in the row I* x with the same coefficients as 

in the massive s-holomorphic row-to-row propagation Pp, in Lemma 6. But by Theorem 10, these 
linear combinations are just the inverse induced rotations applied to tp in the row I*_i, i.e. the 

y~2 



DISCRETE HOLOMORPHICITY AND ISING MODEL OPERATOR FORMALISM 



21 



definitions of the fermions ip on horizontal edges in the row I* j . Again ifi is recovered by the 

application of J. This proves the existence of the extension satisfying the local relations (4.2) and 
(4.3). □ 

4.3. Ising parafermionic observables and low temperature expansions. 

4.3.1. The two-point Ising parafermionic observables. We next consider graphical expansions of cor- 
relation functions of the fermion operators ip(z), tp(z). These are expansions in powers of the 
parameter a = e~ 2/3 , and they are called low temperature expansions because the parameter is small 
when the inverse temperature is large (a \ as f3 /* oo). 

Let a G I* x J be a horizontal edge and z £ (I* x J) U (I x J*) any edge of the rectangle 
I x J. The set of faces I* x J* of the rectangle forms the dual graph, and we denote by 8* = 
{< p,p' > | p,p' € I* x J*, \p — p'\ = 1} the set of dual edges. The low temperature expansions of 
fermion correlation functions will be simply expressible in terms of the following two parafermionic 
observables: 

f\{z) = \ £ ^)e-^) 

/*(*) = ! £ ^e-iW^M, 
-reci(z) 

where the notation is as follows: 

• Ca{z) is the set of collections 7 C £* of dual edges such that the number of edges of 7 
adjacent to any face p G (I* x J*) \ |a + i|, pi 7 ^ is even, and the number of edges adjacent 

to a + \\ and p^ is odd, where p^ is one of the faces next to z. The set Ca{z) is defined 

similarly, but the exceptional odd parities are now at a — i| and at p^ one of the faces next 
to z. We visualize 7 as in Figure 4.1 as a set of loops on the dual graph, together with a 
path from a to z starting upwards/downwards from a, by including two "half-edges": from a 

to a ± i?j and from p^ to z. 

• For 7 G C a (z) we let £(7) = (7] + 1 denote the total length of the loops and the path, 
where (7! is the cardinality of 7 C £* and the additional one is included to account for the 
the two half-edges. 

• The number W(7 : a — > z) is the cumulative angle of turns along a path in 7 from a to z. 
The path is not necessarily unique, but if it is chosen in such a way that no edge is used 
twice and no self-crossings are made, then one can show that the winding is well defined 
modulo 47r and thus the factor e _ 2 ^(i-^z) - g we jj jggjjgj [HoSmlOb]. 

• Z is given by Z = ^ weC a' u ', where C is the set of collections u C £* of dual edges such 
that the number of edges of u adjacent to any face p 6 I* x J* is even. We visualize ui as a 
collection of loops. The expression for Z is the low-temperature expansion of the partition 
function, and it is easy to see that Z = Z + x const., where the constant is e / 3x l( I * xJ ) u ( IxJ *)l. 

The Ising parafermionic observables are s-holomorphic and satisfy the Riemann boundary condi- 
tions, with a discrete singularity at z = a. To give a more precise statement, we first define a notion 
of discrete residue. 

Definition 20. Let a be a horizontal edge. For a function z 1— > f(z) that is (massive) s-holomorphic 
for z 7^ a in a domain containing the faces a± |, the discrete residue of f at a is Res a (/) = 
_L ^front^ _ yback^ a ^ w h ere / front (a) is such that if / is extended to a by this value, then / 
becomes (massive) s-holomorphic on the face a + 5, and / back (a) is such that if / is extended to a 
by this value, then / becomes (massive) s-holomorphic on the face a — ^. 
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a 



FIGURE 4.1. A configuration in Ca(z). The winding in the picture is W(7 : a — > z) = —2ir. 

Proposition 21 ([HonlOa]). Let a G I* x J. If a is not on the boundary, a G I* x (J \ <9J) 7 then 

a I 

i/te Jsmj parafermionic observables f a and fa are functions defined on edges z ^ a such that 

• z i — y fl(z) and z i— > are massive s-holomorphic 

A. I — I I 

• /a and /o satisfy RBVP: for z a boundary edge of the rectangle f a (z) G Mr cw 2 and fa(z) G 

• i/ie discrete residue of fa at a is ^ and the discrete residue of fa at a is 

If a is on the bottom boundary, a G Iq, i/ien /a is ^ero and /„ is a function defined on edges z ^ a 
such that 

• z i — y fl(z) becomes s-holomorphic in the whole domain with the definition /2(a) = 1 

• fa and satisfies RBVP: for z a boundary edge of the rectangle f a (z) G Mr cw 2 . 
If a is on the top boundary , a Gl* N , similar statements hold. 

The parafermionic observables can be defined similarly in any square lattice domain [HoSmlOb]. 
At the critical point, (3 = j3 c , one can treat scaling limits as follows. Take the domains to be 
subgraphs Q of the square lattice 57? with small mesh 6, approximating a given continuous domain 

as 5 \ 0. The analogue of the above Proposition holds. The convergence of the parafermionic 
observables as 5 \ can be controlled [HoSmlOb, HonlOa]: the functions fa and fa divided by 5 
converge uniformly on compact subsets of J7\{a} to the unique holomorphic function with Riemann 
boundary values and the appropriate residue. By Theorems 22 and 23 below, we can deduce from 
this also the convergence in the scaling limit of the renormalized fermion correlation functions. 

4.3.2. Fermion operator two-point correlation functions. 

Theorem 22. We have 

M*Ma)) IXJ = -fl(z) + xfi(z) 
(VW(a)} IxJ = fl(z)+ifHz) 

= - J[a)-ijHa) 
<^(^(a)) IxJ = -Jiz)-iJ{i). 

Proof. Note that because of the relations of Theorem 19 and the fact that f a are massive s- 
holomorphic, it suffices to prove the statements when z is a horizontal edge. Denote z = x + iy 
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and a = x' + iy' . Suppose for simplicity first that y > y' . Consider the numerator of the second 
correlation function, 

{e (+) \v N -^ x vy-y'^ x ,vy'\e {+) ) 

Expand the matrix product in the basis (e a ). Note that for any given a G {il} 1 the matrix elements 
(Vv) Tcr an d (V'aOro- are non-zero only if r is obtained from a by flipping the spins on the left of x' 
or x. The expansion is 

(e i+) \v N ~y^ x vy~y'^vy'\e M ) 

= COnSt. X ^ V( + ), a (N-l) V a (N-l) , a (N-2) V a (N-2) , a (N-3) ■ ■ ■ 
" " " K-d/+l),T(w) (^riy)^) V a (y)^(y-1) ■ ■ ■ 

■ ■ ■ V a(y'+ 1 ),r(y') (Vv) T („') j(7 („') ^(s'),^'-!) • • • 

where the sum is over indices G {±1} and the flipped spin configurations are 

for j > x (y') j ^ ^ or 3 ^ x ' 

for j < x. j |— of ^ for j < x'. 

For the matrix elements of V use the formula 



j 




= g f E fce i*P fc _iP fc+ i xe /3E je i^'P J x e 2 E fce i*o- fc _i^ fc+ i 

= const, x a fc 2 fc +3 1 J ^ Pjt a 2 fc +2 



where a = e _2/3 . In most rows, we can combine the factors from the matrix elements of two (V h ) 5 
to just one factor. The sum essentially amounts to summing over spin configurations in the entire 
box, except from the peculiarity that in rows y and y' we have two configurations related to each 
other by flipping the spins on the left of x or x'. Thus the terms in the sum correspond to contours 
7 G Ca(z) U Ca{z) by the rule that a dual edge is in 7 if it separates two spins of opposite value: in 
rows y and y' the two flipped configurations amount for half-edges arriving to the points z = x + iy 
and a = x' + iy' . The half edge in row y has two possible directions. The half-edge is either from 
x + iy to the face x + i(y + i) above (resp. the face x + i(y — I) below) if a^ y K = a^ y \ and 

/ (resp. cr^i / cr^i and r^i = T^ y \) and in this case we set 77 = +1 (resp. 77 = — 1). 

x '2 3:1 2 :E +2 21 2 ^+2 21 2 

Similarly we set rf = +1 or rf = — 1 if the half edge in row y' is from x' + iy' to the face above 
or below, respectively, i.e. if a^ v \ = a^ y \ or a^ y \ / &^ y ,\, respectively. The matrix elements 

X +2 X j x +2 ^ 2 

of all V together produce a factor a L< - 7 ^ times a constant. The matrix element of ip x produces the 
complex factor i(— l)^^ nI y x ^ , where #(7 n 1^ x ) is the number of edges of the contour 7 on row 
y on the right of x and A = e l7F//4 . Similarly the matrix element of ip x i produces the complex factor 

(— 1) y' X~ v . We now write the result of the expansion in terms of sum over contours, 

{e {+) \v N ~y^ x vy-y'^ x ,vy'\ G{+) ) 

= const, x i £ a ^(7) ( _i)#W)+#(^-') A W. 
7 ecl( 2 )uci(2) 

Combinatorial considerations of the topological possibilities for the curve in 7 from a to z show that 
(_l)#(T nI ^)+#(T nI ^ )\v-rf = _ ie -iW( 7 :a^) ) where W ( 7 . Q _^ ^ ig the win( jing of the path as 
in the definition of the parafermionic observable (note a difference to the case y < y': we would 
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have (_ 1 )#(T nI ^)+#(7m^ ) xv -rf = ie -^W( r .a^z) instea d). Thus we write our final expression for 
the numerator of the second correlation function, 

(e i+) \v N -y?p x vy-y'ijj x ,vy'\e i+) ) 

= const, x « L(7) e-2 w (^). 

^Cl(z)UCi(z) 

The denominator is (e^\V N \e^) = const, x Z with the same multiplicative constant (here Z is 
as in Section 4.3.1), so we get the expression 

{^{a)) = fl{z) + xfi(z). 

In the case y < y' , before we do the expansion of the matrix product, we must anticommute ip(z) 
to the right of tp(a), which gives an overall sign difference. This is nevertheless cancelled in the 
end result by another opposite sign resulting from the combinatorial considerations of topological 
possibilities for the curve 7. 

For the first correlation function, a similar consideration gives when y > y' , 

(e i+) \v N -y^ x vy-y'^ x/ vy'\e {+) ) 

= const, x (-1) £ ^(7) ( _ 1) #(7m-) + #(7ni-') A ^_ 

■yed(z)uci(z) 

In this case we have (_l)#T n V + #T nI / = ^-iwf^^ leading to 

<lK*W(a)>= -fl{z) + xfi(z). 

□ 

4.4. Pfaffian formulas for multi-point fermion correlation functions. The multi-point cor- 
relation functions of the fermions can be written in terms of two-point correlation functions. Recall 
the abbreviated notation of Section 4 for fermion correlation functions — in particular each i/j® in 
the statement below can be either ip or ■0. 

Theorem 23. We have 



^ l \z 1 )---^ n \z n )) =Pf 
/ IxJ 



IxJ 



Proof. We use the polarization of Lemma 15, which works for all (3 > except possibly isolated 
values, and since both sides of the asserted equation are analytic as functions of ft, the statement 
will be proven for all (3. By the aforementioned lemma, the state v vac = e( + ) is a vacuum of the 

Fock space S + = f\W^' N , and the mapping 



e (+) V e (+) 



defines the dual vacuum w* ac G (A ^''^j ■ The denominator in the definition of correlation 
functions in Section 4 is the same as the denominator in the above formula for the dual vacuum, 
(e^lV^le^)) = eJ + ^V N e^ + y The correlation functions thus read (ip^^zi) ■ ■ ■ ^ n \ z n))i x j = 

(v* ac , ip^\zi) ■ ■ ■ ip( n \z n )v vax: ) . Finally note that ip^(zi) G W for all i = 1, 2, . . . , n, so the state- 
ment follows from the fermionic Wick's formula, Lemma 12, applied to the polarization W = 
W^ );t of Lemma 15. □ 
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4.5. Multipoint Ising parafermionic observables. Let us now define multipoint parafermionic 
observables introduced in [HonlOa]. Let be a square grid domain, with dual £1* consisting of 
the faces. Denote the set of edges of $7 by £ and the set of dual edges by £*. Let Z\, . . . , Z2 m be 
(midpoints of) edges, and for each Zj, let Oj be a choice of orientation of the corresponding dual 
edge e* (i.e. Oj G {±1} if e* is horizontal and Oj G {±«} if e* is vertical), and let Ej G C be choices 
of square roots of the orientations, ej = oj. 

We define the multipoint observable f e (z\, . . . , Z2 m ) by 

r(z 1 ,...,z 2m )= £ assign ( 7 ) ]J f* e -*W0*>, 

where 

• C| 1; _ Z2m is the set of 7 C £* consisting of the (dual) half edges < Zj, zj + ^ > and of (dual) 
edges of £ * such that each vertex p G $7* belongs to an even number of edges /half edges of 
7: in other words a configuration 7 contains loops and m paths tt±, . . . , ir m linking pairwise 
the Zj's. By £(7) we mean the number of edges of £* in 7 plus m, with the additional m 
accounting for the 2m half edges. 

• The product is over the m paths 717, . . . , ir m , where each ttj is oriented from z Sj to z&. where 
Sj < dj (i.e. we orient the paths from smaller to greater indices). 

• sign (7) = ^_^^# crossm s s ) w h ere ^crossings is the number of crossings of the pair partition 
{{sj, dj} : j G {1, . . . , m}} of {1, . . . , 2m} induced by the paths m, . . . , 7r m (irj from z Sj to 
z,£.), i.e. the number of 4-tuples Sj < dj < < d^. 

• It can be checked [HonlOa] that if there are ambiguities in the choices of paths 717, . . . , n n , the 
weight of a configuration 7 is independent of the way that they are resolved, provided that 
wherever there is an ambiguity each path turns left or right (going straight is forbidden) . 

The observables f e (27, . . . , Z2m) can be used to compute the scaling limit of the energy density 
correlations, as well as boundary spin correlations with free boundary conditions (see [HonlOa]). 
The key property that allows one to study the observable at criticality is its s-holomorphicity: 

Proposition 24 ([HonlOa]). Let o±, . . . , o 2m -\ G C be orientations of edges e\, . . . , &2m-\ an d ^ 
£1, ■ ■ ■ , £2m-i £ C be such that e\ = o\, . . . , elm-i = °2m-i- Let z\, . . . , Z2 m -i be the midpoints 
of el, . .. ,e 2 ' m _ 1 . For any midpoint of edge Z2 m , let 02m and 02m be its two possible orientations, 
let £2m and £2 m be such that e\ m = 02m and e\ m := —02m, and let e := (ei, . . . ,e 2m ) and e := 
(ei,... 

Then we have that g (z2 m ) '■= -^f £ (zi, ■ ■ ■ , Z2m)+^/^ (zi, ■ ■ ■ , z 2 m) is independent of the choice 

_1 

of £2m,£2m and at criticality Z2 m ^ 9 (z2m) is s-holomorphic on Q \ {z±, . . . , Z2 m -i}> with || r cw 2 
boundary conditions. 

As for the fermion operator two point correlation functions and two point parafermionic observ- 
ables, it is true that the fermion operator multipoint correlation functions are expressible as linear 
combinations of the multipoint parafermionic observables and vice versa. 

Theorem 25. Define iffi(z) = \{^{z) — ip(z)) and ip^(z) = ^{^p{z) + 4>{z)). Then we have 

U hm {z2m) ■■■^ 1 (zi)) + =f( Zl ,..., Z 2m ), 
\ I Ix J 

where the arrows £j€ {til} and the square roots of directions Ej = ^Jof are chosen as follows 



Ui = A if tj=t 

\ej = A- 3 if tj=i 
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Proof. Suppose for simplicity that Zj = Xj + iyj with y± < yi < ■ ■ ■ < yi m . Form a low temperature 
expansion of the fermion operator correlation function as in the proof of Theorem 22. Consider 
any fixed j = 1,2, ... ,2m. Adding up the low temperature expansions of the two terms in the 
definition of tp$ j (zj) we get that the total weight for the configurations where the half edge from Zj 
to Zj + \r\j\ is used (m G {±1}) is used is in the two cases Xj = t an d Xj = X respectively proportional to 
^(A - ^ — iA^') = \~ 1 5 Vj:+ \ and ^(iA^ + A - ^ ) = A 3 (5 %i _i, where the Kronecker deltas in particular 
ensure that only the contributions of the contours Cf _ survive, as in the definition of the 
corresponding parafermionic multipoint observable. In the low temperature expansion, contours 7 
always come have a factor a L W in their weight due to the product of matrix elements of V h and 
V v , and for any surviving contour 7 G C Z1 Z2n the remaining phase factor coming from the matrix 

elements of (ip\ 3 ^ ec l ua ls sgn(pairing(7)) x \Y^ =l < j s -. □ 

As a direct consequence of Theorems 23 and 25, we get a PfafHan formula for the multi-point 
parafermionic observables. 

Corollary 26. Let e\, . . . , e\ n be dual edges with orientations 01, . . . , 02 n and let e±, . . . , Em be such 
that e\ = 01, . . . , e\ n = 02 n - Then we have that 

f (z u z 2n ) = Pfaff (f^^ ( Zj , z k ) l# k ) 

This formula was proved in [HonlOa] in the critical case /5 = (3 C for domains of arbitrary shape, by 
verifying that the s-holomorphic function in Proposition 24 satisfies a discrete Riemann boundary 
value problem with singularities, which uniquely characterizes the parafermionic observable. The 
special case which gives the Ising model boundary spin correlation functions with free boundary 
conditions was proven by direct combinatorial methods in [GBK78, KLM12] for very general classes 
of planar graphs. Our approach works at any /3, but for the Ising model on square lattice only, and 
for domains of general shape some minor technical modifications are needed in the proof: the domain 
should be thought of as a subgraph of a large rectangle IxJ, and for every row the transfer matrix 
should be replaced by a composition of V and a projection which enforces plus boundary conditions 
outside the domain. Nevertheless, we believe that our approach in conceptually the clearest, as the 
PfafHan appears simply because of the fermionic Wick's formula (Lemma 12). This illustrates an 
advantage of the operator formalism, some algebraic structures underlying the Ising model are more 
evident and can be better exploited. 

4.6. Correlation functions of the fermion and spin operators. It is also possible to consider 
correlation functions of fermion operators and spin operators simultaneously. It turns out that as 
functions of the fermion operator positions, these become branches of multivalued observables. For 
example, when a G Iq is on the bottom side of the rectangle and w±, . . . , w n G I x J, 

{ip{z)i){a)a{wi) ■ ■ ■ o-(w n ))i xJ _ (e (+) \V N il>(z)4>{a)d-{wi) ■ ■ ■ a(w n )\e (+) ) 
(a( Wl ) • • • 6-{w n )}f xJ (e {+) \V N a(w 1 )---a(w n )\e {+) ) 

becomes a (massive) s-holomorphic function in the complement of the branch cuts starting from 
each Wj G I x J to the right boundary of the rectangle. The function can be extended to the branch 
cut in two ways, one of which satisfies (massive) s-holomorphicity conditions on the faces below 
the cut and another which satisfies them on the faces above the cut — the two definitions differ 
by a sign, indicating a square root type monodromy of the function at the locations Wj of the spin 
insertions. A low temperature expansion like in the proofs of Theorems 22 and 25 shows that this 
function is a branch of the parafermionic spinor observable of [Chlzll, CHI12], where the observable 
is properly defined on a double covering of the punctured lattice domain in order to obtain a well 
defined s-holomorphic function. 
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5. Operators on Cauchy data spaces 

As explained above, the Ising transfer matrix can be constructed directly in terms of the s- 
holomorphic propagator (Sections 1.3.1 and 3.3). We now discuss s-holomorphic approaches to the 
data carried by the transfer matrix quantum states. In Sections 4.3 and 4.5, we learned the following: 

• The correlation functions of the fermion operators can be expressed as linear combinations 
of parafermionic observables. 

• The parafermionic observables can be characterized in s-holomorphic terms: they are the 
unique s-holomorphic functions with Riemann boundary values and prescribed singularities. 

In this section, we present an s-holomorphic construction inspired by transfer matrix states. A 
quantum state Q 6 5 living on a row contains all the information about the geometry of the 
domain and the operator insertions below Likewise, we construct discrete Riemann Poincare- 
Steklov (RPS) operators living on a row (more generally, any crosscut of the domain), which act 
on Cauchy data spaces. These RPS operators together with the vectors on which they act contain 
all the information about the geometry of the domain and operator insertions. These operators can 
be written as convolution operators with parafermionic observables, which are fermion correlations 
and hence can be directly represented from the quantum states. They also can be propagated using 
explicit convolution operators. 

A great advantage of the discrete RPS operators is that they have nice scaling limits, as singular 
integral operators, and that they work in arbitrary planar geometries. 

5.1. Discrete RPS operators. In this Section, we define the Riemann Poincare-Steklov operators. 

Let SI be a square grid domain, let b C d£ be a collection of boundary edges and and be the 

_ i 

space of functions / : b — > C such that / || r ccw on b. Let be the space of functions / : b — > C 
_ i 

such that / || Tew 2 on b. 

First, we state a key lemma, which guarantees the uniqueness of solutions to Riemann boundary 
value problems. 

Lemma 27. Let Q be a square grid domain with edges £. If h : £ — > C is an s-holomorphic function 
_ i 

with h || t cw 2 on d£, then h = 0. 

Proof. The proof of this lemma is given in [HonlOa, Corollary 29] (where the notion of s-holomorphicity 
comes with a phase change of e l7r / 4 compared to the present paper) . The idea is to show that for any 
s-holomorphic function g : £ — > C with boundary values u + v, where u G TZq and v Gl^, we have 
that XLeb \ v ( z )\ — Szefa \ u ( z )\ ( tne P ro °f °f this inequality relies on the definition of a discrete 
analogue of 9m J g 2 ). In our case u = 0, and hence v = as well. □ 

Lemma 28. For any u G TZq, there exists a unique such that u + v has an s-holomorphic 

_ i 

extension h : £ — > C satisfying h \\ r cw 2 on dfl \ b. 

Proof. For any v G Iq, there exists at most one u G TZq such that u + v has an s-holomorphic 

_ i 

extension h to SI — > C with || r cw 2 on d£ \ b: if we suppose there are two extensions, their difference 

_ i 

will satisfy the boundary condition || r cw 2 on d£ and hence be by Lemma 27. By a dimensionality 
argument, there exists exactly one such u and the mapping v i— > u is an invertible linear map. □ 

Definition 29. We define the RPS operator Uq : TZq — > Iq as the mapping u >->■ v defined by 
Lemma 30 (which is an isomorphism by the proof). 

Lemma 30. With the notation of Lemma 28, we have that 

(5.1) v(x) = Yl u(y)f n (y,x) Vx G d£ 

yeb\{x} 
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(5.2) h{x) = Y J u{y)fa{y,z) Vxef. 

Remark. The convolution formula, Equation (5.2), is the key to pass to the scaling limit: the kernel 
fn converges to a kernel of a continuous singular integral operator (see [HoKyll, Section 13]). 

_ i 

Proof. Let us first notice that (5.2) implies (5.1): if x € d£, we have u (x) fn (x, x) £ r CC w and 
i 

u (y) fn (y, x) £ t cw 2 (y) for y £ d£ \ {x}, and hence the projection of (5.2) on I is indeed (5.1). 

To prove (5.2), notice that the right-hand side is an s-holomorphic function, with the right bound- 
ary conditions (i.e. the same as h in Lemma 28). The difference of both sides must then be 0, by 
Lemma 27. □ 

In rectangular boxes, the RPS operator can be written simply in terms of the s-holomorphic 
propagation 

Lemma 31. Let £1 be a rectangular box I x {0, ...,N}, let b = I X {0} be the bottom side and 
let P : (M 2 ) 1 — > (M 2 ) 1 be the s-holomorphic propagation as defined in Section 2-4- For N > 0, 
decompose the s-holomorphic propagation P N into four |I*| x |I*| blocks 



, pN pN 
r — 1 pN pN 



corresponding to the decomposition (M 2 ) 1 = M 1 © iM 1 into real and imaginary parts. Then have 

Uh = - {Pqq) Pm- 

Proof. Let u 6 TZq (i.e. purely real in this case) and v G I Q be defined by v = UqU. By definition 
of Uq, we have that 

for some purely real function w : I x {n} — > M. Hence, we get that PqshU + P^v = 0. Since we know 
that for any u G TZ n , there exists a unique v satisfying this equation (as is an isomorphism), we 
get that v = - {P^)' 1 P^u. □ 

5.2. RPS pairings. In this subsection, we explain how to pair together s-holomorphic data coming 
from two adjacent domains with disjoint interiors, using RPS operators and Ising parafermionic 
observables. A related discussion about discrete kernel gluings (in a different framework, without 
boundaries) can be found in [Dublla]. In our framework, the gluing operation arises as an analogue 
of pairing of transfer matrix states. 

Let us define the setup of this subsection. Let ^1,^2 be two adjacent square grid domains with 
disjoint interiors, with edges £i,<?2, let ft := ^1 U ^2 an d assume that b := d£\ Pi 882 is connected. 
Let U\ := Uq and U2 '■= Uq 2 be the RPS operators defined in the previous subsection and set 
Kj := K b n . and X) := 1^. for j = 1, 2. We have TL X = X 2 and H 2 =T X . 

Lemma 32. We have that (Id — U-JJ2) '■ TZ<2 7^2 an d (Id — U2U1) : 1Z± — > 1Z± are isomorphisms. 

Proof. The injectivity (and hence the bijectivity) of these operators follows from the fact that if a 

function u £ IZ2 is a fixed point of U1U2, then u + U2U admits an s-holomorphic extension to f2 with 

_ 1 

boundary condition || r cw 2 on d£, and is hence by Lemma 27. □ 
A useful corollary of the previous lemma is the following fixed point result: 
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Corollary 33. Let h\ G X\ and h 2 G Z 2 . Then there exists a unique function f : £ — > C such 

that for j = 1,2, the function f — hj has an s-holomorphic extension to £j with boundary conditions 
_ 1 

|| r cw 2 on We have that f = u± + U2, where u\ and u 2 are given by 

ui = (Id - C/2^1)" 1 (U 2 hi + h 2 ) , 
u 2 = (Id-C/i?7 2 )- 1 (?7 1 /i 2 + /i 1 ). 

Proof. Suppose first that there exists an / such that the functions / — hj have s-holomorphic exten- 

_ 1 

sions to £j with boundary conditions || r cw 2 on d£j. Set fj := f — /ij and write /j = Uj + fj, where 
G 7£j and G Xj . We have that f = ui + u 2 and 

hi = f 2 + /12 = ^2^2 + h 2 = U 2 (vi + hi) + h 2 = U 2 U\ui + U 2 h\ + /i 2 
^2 = v\ + h\ = U\ui + hi = Ui (v 2 + h 2 ) + hi = UiU 2 u 2 + f/i/i 2 + hi, 

which gives that (Id — U 2 Ui) ui = U 2 hi + h 2 and (Id — UiU 2 ) u 2 = Uih 2 + hi. By Lemma 32, we 
obtain the asserted formulas for ui and u 2 , proving the uniqueness of /. For any (hi, h 2 ) we have 
seen that there is at most one and hence exactly one solution (ui,u 2 ) of the equations 



□ 



' Ui 


-id " 




Ul 




' -hi ' 


-id 


U 2 




u 2 




-h 2 



showing the existence of / with the desired properties. 



As illustrated in the next subsection, Corollary 33 has the following consequence: the value (on 
b) of an s-holomorphic observable h : Q — > C with Riemann boundary conditions and prescribed 
singularities in fii, Q, 2 can be recovered from (hi, Ui) and (h 2 , U 2 ). In other words, these pairs carry 
all the relevant information about fii, f2 2 that is needed to compute s-holomorphic correlations. 
More precisely, and as will be illustrated in the next subsection: Ui,U 2 encode the geometry of the 

domain and hi,h 2 encode the singularities (in practice, they are the restriction to b of functions 

_ 1 

with singularities in f2i,f2 2 and r cw 2 boundary conditions on d£i,d£ 2 ). 

Once the values of an s-holomorphic observable h on b are known, one can compute the values of 
h — hi on £ 1 and the values of h — h 2 on £ 2 , using the convolution formula (5.2) in Lemma 30. 

5.3. Fermion correlation and fixed point problems Cauchy data. The fermion correla- 
tor/parafermionic observable fits naturally in the framework of the previous subsection. A first 
consequence is the following. 

Proposition 34. With the notation of Section 5.2, set Q := (Id — UiU^ 1 . For any x £ d£i \ b, 

we have 

fn(x,-) r =(ld + U 2 )Qf Ql (x,-) 



Proof. Set / := fa (x, •) , fi := f Ql (x, •) 



Write / = m +u 2 where Uj £ Hj. Applying Corollary 



33 to /, hi := fi and h 2 = 0, we get u 2 = Qfi- Since fa (x, •) is s-holomorphic with r cw 2 on d£ 2 \ b, 
we have that / = u 2 + U 2 u 2 . □ 

Extending fa(x,-) to £ 2 gives in particular the following nice formula: 

Theorem 35. For any x G d£i \ b and y G £ 2 , we have 

T I 

(5-3) fa(x,y) = fa 2 (-,y) Q fa 1 (x,-)\ 

b lb 
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Proof. Set / := fa (x, ■) 



and write / = 112 + V2, with U2 G TLi and U2 G I2. By Lemma 30, we 



have that fa(x,y) = J2 zeb u 2 (z) fa 2 (z,y). By Proposition 34, u 2 = Qfa, (x, ■) 
follows. 



and the result 
□ 



The formula is an analogue of a natural pairing in transfer matrix formalism: for example when 
$7 = I x {0, 1, . . . , N}, and x G Io, x' + iN G Ijy, the correlation function (ip (x) if) (x' + iN)} can 
be obtained by propagating the state ip x i with V k to the fc:th row and pairing it (with the inner 
product in S) with the state ^J,f propagated downwards from row iV to row k: Equation 5.3 is the 
analogue of 

When x and y are both in E\ , we can pair the states associated with x and y (this is not possible 
with transfer matrix): 

Proposition 36. For x G d£\ \ b and y G S\, we have 

fa (x,y) = fa, (x,y) + (fa, (;y) \ b f (Id - U^y 1 U 2 (fa, (x, •) | J . 

Proof. Set f = fa (x, •) and write f = u\+ v\, with u\ G 1Z\ and v\ G Xi. By Lemma 30, we have 

fa 

that 

fa (x, y) - fa, (x, y) = ^2,ui (z) fa, (z, y) 

zeb 

As in the proof of Proposition 34, with Corollary 33 applied to h\ = fa r (x, ■) 
that u\ = (Id — UiU2)~ 1 U2 (fa, (x, •) ^ and the result follows. 



and h 2 = we get 



□ 
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